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Abstract 



, The basic equations describing a Newtonian universe with uniform pressure 

are reexamined. We argue that in the semi-classical formulation adopted in 
PLh' the literature the continuity equation has a misleading pressure gradient term. 

Q ■ When this term is removed, the resulting equations give the same homoge- 

neous background solutions with pressure and the same evolution equation 
\ for the density contrast as are obtained using the full relativistic approach. 
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I. INTRODUCTION 



Seventeen years after Einstein's seminal paper on relativistic cosmology, Milne and Mc- 
Crea [|I|] obtained Newtonian analogs to the expanding dust-filled FRW models of the uni- 
verse, thereby initiating the so-called Newtonian cosmology. Apart from some conceptual 
ambiguities present in these pressureless models, there is no doubt that such analogs are 
relevant both from a pedagogical and a methodological point of view, since they describe, 
in a rather simple way, several features of the presently observed universe. A natural gener- 
alization of these models should include pressure effects, as an attempt to construct analogs 
to the full class of FRW models. Unfortunately, such an extension cannot be analytically 
implemented, at least when the approach is based uniquely on Newtonian concepts (absolute 
time, Euclidean space and ordinary gravity). The basic difficulty arises because a uniform 
pressure, in the Newtonian framework, does not play any dynamical or gravitational role, 
at the level of the continuity, Euler and Poisson equations. A "neo-Newtonian" description, 
however, has been derived long ago by McCrea himself, and in a more definite form 
(without any concepts from general relativity), by Harrison ||^. The later approach is based 
on the following set of equations: 

"Continuity" equation: 

| + V;.(p + ^)n = 0, (1) 

"Euler" equation: 

— +U.VrU = -Vr^-ip + ^)-WrP, (2) 

at & 

"Poisson" equation: 

V^^ = 47rG(p + 3^) , (3) 

where c is the velocity of light and the quantities p, p, u and tl) denote, respectively, the 
density, pressure, field velocity and generalized gravitational potential of the cosmic fiuid. 
As discussed in detail by Harrison 0, all corrections of order are based uniquely on 
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arguments coming from special relativity. The price to pay is that the standard Eulerian 
equations are recovered only in the limit c oo. However, the main pedagogical aspect 
contained in the original equations for a dust-filled universe has been preserved, viz., all 
the fundamental Newtonian notions like absolute time, Euclidean space and gravity are also 
part of the enlarged theoretical context represented by the above modified equations. 

II. A NEW "CONTINUITY" EQUATION 

As is widely known, the above set of equations admits a homogeneous and isotropic 
solution, i.e., p = p{t) and p = p{t). In this case, the fluid velocity is (a dot means partial 
time derivative) 

u = - f, (4) 
a 

where the evolution of the scale factor a{t) is governed by the Friedmann equations, namely: 

SttG kc^ . . 

with the continuity equation (|1|) reducing to the following form 



9p .^a p 
— + 3-(p+ - 
ot a c 



+ 3-(p+^)=0. (7) 



Equations (5)-(7) are exactly the same as obtained from Einstein's field equations, when 
one adopts comoving coordinates to describe a homogeneous and isotropic spacetime filled 
with a perfect simple fluid. Paradoxically, if one studies perturbation theory in the frame- 
work of equations (|lD-®, the evolution equation for the density contrast does not agree 
with the relativistic result in synchronous gauge. The expected result is obtained only for 
a dust-filled universe {p = 0). This problem was first mentioned in a short note written in 
the appendix of the paper by Sachs and Wolfe M. 
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The root of the problem seems to be closely related to equation (|T|). As we shall see, this 
equation corresponds to the correct conservation equation only in the homogeneous case. If 
inhomogeneities are taken into account, as for instance in the theory of small fluctuations, 
the pressure gradient appearing in Eq.(l) generates extra terms, leading to results which 
disagree with those resulting from the relativistic theory. To clarify this point, let us consider 
the energy conservation equation in the form dE + pdV = 0. The energy contained in a 
volume V with radius r ~ a(t) is ~ pc^a^, and since the volume is changing we get 

dp , P . 1 d n 

Using Eq. (|D we have Vr-u = 3d/ a, and then it follows that 

Finally, by expanding the total derivative, = §i + u.Vr, the above equation may be 
rewritten as 



| + V..(,3)H-| 



+ V,.(pm) + ^V,.M = 0. (8) 



This is the correct continuity equation in the modified Newtonian approach to cosmology 
when pressure effects are included. It looks like a balance equation for the matter density. 
The source term (the last term in Eq. (^), is related to the work dW = pdV needed to 
expand the fluid from the volume to + dV. In fact 

1 dW Ana'^da d ^ ^ 

T7-Jr= P~i — = = P^^-'^ ■ 
V dt -^Tra-^dt a 

where we used Eq. in the last equality. This term has to be added to the usual equation of 

fluid dynamics to account for the work related to the local expansion of the fluid. It is worth 

noticing that in the limit of uniform pressure, p = p{t), Eq. (|^) can always be rewritten in a 

form resembling Eq. (0), the form usually adopted in the literature [?] (Note that McCrae 

[?] originally wrote down the equation in the form of §!). However, although formally 

reducing to the same conservation equation for a homogeneous distribution, the perturbed 

equations that follows from each case are quite different. As we shall see, unlike Eq. (|l]), the 

new "continuity" equation (P) is completely consistent with relativistic perturbation theory. 

In particular, this solves the contradiction pointed out by Sachs and Wolfe. 



III. PERTURBATION THEORY 



To study the evolution of small fluctuations in an expanding Newtonian universe, let us 
now consider the standard perturbation ansatz (in this section we take c = 1): 

p = p,(t)[l + 5(f,t)], (9) 
p = Po{t) + 6p{f,t) , (10) 

^ = ^o(f,t)+^(f,t), (11) 

u = Uo + v{f, t) , (12) 

where the quantities carrying the subscript "o" represent the (background) homogeneous 
solution to the unperturbed equations (0), (||) and (H). The quantities 5, 5p, ip and v are 
to be considered small (perturbations) when compared to the their respective background 
quantities. It should be recalled that in terms of the comoving coordinates q we have 
r{t) = a(t)q and Uo = a(t)q. 

Inserting the above expressions into Eqs. (0), (|^) and (§), and linearizing the resulting 
equations, we get to first order in the perturbations 



Po 



di 



— + Uo-VrS 



3^Po5 + 3^5p+{p + p)Vr.v = 0, (13) 



— + {Uo.S/r)v+ -V = -V,V^ - {p + pY^S/ r6p , (14) 

ot a 

V^if = 47rGpo{6 + 3—) . (15) 

Po 

We now change to comoving coordinates q = f/a. Following standard lines [^], we 
transform the partial time derivative of an arbitrary function / at fixed r,^ , to the partial 
time derivative of / at a fixed q, which we write as / . The relation between these two partial 
time derivatives is 

where V = aVr is the gradient with respect to q at fixed time. 
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For the sake of simplicity, let us also assume that the medium satisfies the equation 
of state p = up , where u is a constant. For adiabatic perturbations, this means that 
6p = Vg6p = v1po5, where f ^ = z/ is the sound velocity. With this choice and using Eq. (^6|), 
Eqs. (|T3|)-(|1|) reduce to 

5 + ^PoV.t; = 0, (17) 
a 

'v + -v = --V^- W, (18) 

a a (1 + v)q 

VV = 471^(1 + ^v)poa^5 . (19) 

Finally, by eliminating the peculiar velocity from Eqs. (17) and (18), and using (19), is 
readily seen that 

2 

5 + 2-5-47rG'(l + i^)(l + 3i^)5 = %V2(5. (20) 
a a'' 

This is the differential equation governing the evolution of the density contrast, when we 
describe a mass distribution with uniform pressure using the modified Newtonian equations. 
It coincides with the general relativistic equation (in the synchronous and comoving gauge), 
for any value of the parameter u. In particular, for scales bigger than the Jeans' length, 
when the Laplacian term on the right hand side may be neglected, Eq. (^) give rise to the 
same kind of gravitational instability usually derived in the relativistic treatment (see, for 
instance, Eq.(10.118) in Peebles' book ||^). 

2 

As an illustration of these results, we consider the unperturbed solution, a(t) ~ ^3(1+"), 
which is obtained when the field velocity coincides with the escape velocity (fiat model, in 
the relativistic terminology). In this case, po ~ a^^'-^^'^-*, and for scales bigger than the 
Jeans' length, the growing mode of ( pOD scales as 5+ ~ a^+'^^. It thus follows, from (|19|), 



that the correction to the Newtonian potential is time- independent, regardless of the value 
of u. As a consequence, there is no contribution to the temperature anisotropy of the cosmic 
background radiation due to the integrated Sachs- Wolfe effect. In fact, this effect is given 
by 
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ST 



ff 

n) lisw = J dT(^{T, n(ro - r)) , (21) 



where the integral is taken along the photon path from the last scattering surface to the 
observer [^]. Since the perturbed Newtonian potential is time independent for any value of 
it thus follows that ^ |isw = 0. This generalizes the well known result for a dust-filled 
universe. Note that this result (time independent Newtonian gravitational potential) is in 
agreement with what we expect from a Newtonian analysis. 

For completeness, we also derive the evolution equation for the density contrast, using 
the standard "continuity" equation ([1|). Inserting the perturbation ansatz in Eqs. (|l])^(il) 
and following the same steps as above, we arrive at the following equation 

2 

6 + 2-5 - AnG(l + + 3u)5 + %q.V5 + -uq.V5 = . (22) 

a a a 



Comparing to Eq. (pof), we see that (22) has two additional terms (the last two on the 
LHS), and they have no counterpart in the relativistic theory of cosmological perturbations. 
Naturally, their ultimate origin is traced back to the pressure gradient term present in the 
standard "continuity" equation. As a matter of fact, we cannot eliminate these terms by 
simply rewriting the perturbed equations in the comoving coordinate system (where the 
background velocity is zero), as recently claimed by some authors 0. As readily seen from 
(22), such a trick works only for a dust-filled universe(z/ = 0), when both "continuity" 
equations reduce to the ordinary mass conservation law. 

IV. CONCLUSIONS 

We have shown that the correct relativistic equations for the cosmological background 
and for the density contrast can be obtained from Newtonian cosmology even in the presence 
of nonvanishing pressure. This is achieved by a careful consideration of what a modified 
continuity equation should look like. The equation for the density contrast agrees with 
the corresponding relativistic equation written in the synchronous and comoving gauge. 
Moreover, as a byproduct we get a Newtonian potential which is well defined independent of 



the fluid equation of state. In principle, these results imply that the domain of applicability 
of Newtonian cosmology may be enlarged to analyse some problems of structure formation 
even in the radiation dominated phase. 
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